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Constraining the Molecular Gas in the Environs of a z 
Gamma Ray Burst Host Galaxy 
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ABSTRACT 

GRB 090423 is the most distant spectroscopically-confirmed source observed in the 
Universe. Using observations at 37.5 GHz, we place constraints on molecular gas emis- 
sion in the CO(3-2) line from its host galaxy and immediate environs. The source was 
not detected either in line emission or in the rest-frame 850 /im continuum, yielding 
an upper limit of Sgynui— 9.3 ^Jy and M(H2)< 4.3 x 10^ M© (3cr), applying standard 
conversions. This implies that the host galaxy of GRB 090423 did not possess a large 
reservoir of warm molecular gas but was rather modest either in star formation rate or 
in mass. It suggests that this was not an extreme starburst, and hence that gamma ray 
bursts at high redshift trace relatively modest star formation rates, in keeping with 
the behaviour seen at lower redshifts. We do, however, identify a millimetre emission 
line source in the field of GRB 090423. Plausible interpretations include a CO(l-O) 
emitting galaxy at z = 2.1, C0(2-l) at z = 5.2 and CO(3-2) at z = 8.3. Efforts to 
identify a counterpart for the molecular line emitter and to further characterise this 
source are continuing. 
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1 INTRODUCTION 

' Gamma Ray Bursts (GRBs) are amongst the most useful 

' probes of the early Universe. Their utility is two-fold, both 
in yielding information about their host galaxy and imme- 
diate environment, and as beacons, marking out mass and 

, structure at the highest redshifts (Lamb & Rcichart 2000). 

1 Given their high energies, they can potentially be observed 
out to z = 15 or further with current instruments, but red- 
shift determination for these high redshift bursts is chal- 
lenging, requiring deep integrations with optical and near- 
infrared telescopes in the hours after a source is detected. 
Despite this, the few known z > 5 sources have provided 
valua ble evidence for the pr operties of dust ([Stratta et al] 
l-'IMl1 l and molecular clouds jCampana et al.ll2007l ) at these 
early times and for t he colum n density of neut ral hydrogen 
in th eir environment (|Ruiz-Velasco et al . 2007,,; lTotani et al.l 
l2006l ). 

While GRBs do not require massive host galax- 
ies ijSavaglio. Glazebrook. fc LeBorgnj|2009l : [s"vensson et al.l 
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l2010f ). they do require the presence of significant quan- 
tities of gravitationally collapsed and processed baryonic 
matter in the form of young stars. At high redshifts, this 
condition is believed to be sufficiently rare to make them 
likely markers of overdense regions in the cosmic web and 
hence a measure of the power of small sc ale density fiuctu- 
ations (|Mesinger. Perna. fc Haimanll2005l ). Using such bea- 
cons may be the only way to identify and characterize star 
forming regions during cosmic reionization - the epoch dur- 
ing which the first stars heated and ionised neutral hydrogen 
atoms in the intergalactic medium. 

Gamma Ray Burst (GRB) 090423 at z = 8.23+0 07 
is the most distant known object with a well-defined , 
spectroscopically-confi rmed redshift (|Tanvir et al.l l2009l : 
ISalvaterra et al.ll2009l ). Occurring just 600 Myr after the Big 
Bang, it likely probes the tail end of the epoch of reioniza- 
tion. Its afterglow and host galaxy provide an unprecedented 
opportunity to probe the properties of a star forming galaxy 
at these early times, and it is likely to be the archetype 
for future observations on sim ilarly distant GRBs (see the 
Decadal Survey Whitepaper bv lMcQuinn et al]|2009l . for dis- 
cussion). 

Crucially, the properties of GRB 090423 are consis- 
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tent with the bulk of the GRB population jTanvir et alj 
lacking the long duration and extreme energy out- 
put that might be ex pected from a pop III progenitor 
IIMeszaros &: Reej|2010l ). As such, this burst must mark out 
a region in which star formation is already well under way, 
and in which the interstellar and intergalactic medium has 
already undergone processing. For this to have occurred at 
such high redshifts requires this region to have collapsed 
early in the history of the Universe, implying that it is highly 
overdense relative to the volume average at early times. By 
z = 8, it is possible for dark matter haloes as massive as 
a few lO'^^ Mp to h ave collapsed, albeit in small numbers 
||Mo fc WhitJl2002l ). While GRB 090423 is not coincident 
with a known AGN, it could plausibly lie either in a mas- 
sive halo, or in the mass overdensity associated with such a 
halo. 

Interpretation of high redshift bursts, particularly in 
the context of evolutionary trends such as the cosmic star 
formation rate history and inference s made from it (e.g . 
iKistler et al]|2009l : I Wang fc Daill2009l : IWvithe et aLll2010l ). 
requires an understanding of the burst environment and 
how it too changes with redshift. At low redshift, GRBs 
are not unbiased tracers of star formation, but are seen 
to occur preferentially in low mass, metal poor galax- 
ies ([Savaglio, Glazebrook. fc LeBorgnell2009l : ISvensson et al.l 
I2OI0I ). Semi-analvtic models, applied to large simulations 
such as the Millennium Simulation, suggest that GRB host 
galaxy halo mass is largely indepe ndent of redshift to 2: > 6 
l|Chisari. Tissera. fc Pellizza 2010l). but ob servational evi- 
dence for this is sparse. iKistler et M l|2009l ') have suggested 
that the ratio between gamma ray burst rate and star for- 
mation density scales as (1 -I- 2)^'^, possibly due to evolu- 
tion in the mean metallicity. If host galaxy stellar masses, 
star formation rates or dark matter halo masses (and hence 
clustering) of GRB hosts evolve with redshift, then the in- 
terpretation of burst rate as a proxy for star formation rate 
history will need to be revisited and refined. 

GRB 090423 provides a case study at high redshifts, po- 
tentially probing significant evolution in environment from 
2; < 1 burst hosts. Given that the GRB appears to be the 
result of a second generation or later starburst at such early 
times, does it occupy a more massive galaxy than typical of 
lower redshift bursts, a vigourous starburst, a mature sys- 
tem that has exhausted its gas supply or a young one with a 
large reservoir of processed molecular gas available for later 
star formation? 

One test of the latter hypothesis has become possible 
due to the broad bandwidth and high sensitivity of the cur- 
rent generation of millimetre telescopes. At z = 8, the rest- 
frame infrared CO(3-2) emission line (a good probe of molec- 
ular gas mass at low redshift) is shifted into the observed 
7 mm band. The 4 GHz bandwidth available using the new 
CABB correlator at the Australia Telescope Compact Array 
(ATCA) enables a volume with line-of-sight extent Az = 0.9 
to be observed in a single observation, easily encompassing 
the redshift uncertainty from optical/near-infrared observa- 
tions. Given the high redshift of this and similar sources, 
carbon monoxide and other far-infrared emission lines (no- 
tably [CII] at 1898 GHz) may well prove the most accessible 
diagnostics of the precise redshift of the GRB host and any 
galaxies sharing a large scale structure surrounding it. 

In this paper, we describe an investigation to study the 



CO (3-2) emission line properties of the GRB host galaxy 
and neighbouring galaxies in its immediate environs. We 
structure the paper as follows: In section [5] we describe our 
observations, targeted on the host galaxy of GRB 090423 
and searching for CO (3-2) emission. In section [3] we present 
limits obtained on the properties of the gamma ray burst 
host. In sections 14.11 and 14.21 we present the properties of 
additional sources detected in our data, first in continuum 
and then in molecular line emission. Finally in sections[5]we 
discuss the interpretation and implications of our results, 
before presenting our conclusions in section |6l 

All optical and near-infr ared magnitud es in this paper 
are quoted in the AB system (|Oke fc Gunnlll98ay We adopt 
a ACDM cosmology with (^a, ^m, h)=(0.7, 0.3, 0.7). 



2 OBSERVATIONS 

Observations were carried out over the period 2010 Mar 22- 
26th, using the 7mm-band receivers at the Australia Tele- 
scope Compact Arraj0. The 4 GHz bandwidth CABB Cor- 
relator was used, with a small overlap between the two IFs, 
allowing continuous frequency coverage in the range 36.075- 
39.723 GHz, at 1 MHz resolution. A single pointing was ob- 
served at 09''55'°33!29 -M8°08'50'.'80, placing the GRB Host 
Galaxy ten arcseconds north of the phase centre. The pri- 
mary beam FWHM of the ATCA at 38 GHz is 75". 

GRB 090423 is well north of the optimal declination for 
the ATCA. As a result, observations were split into 5 tracks 
and the source was observed at elevations >30° for 4.9hrs 
per track. To ensure good uv-coverage and sensitivity at mil- 
limetre wavelengths, observations were taken in the compact 
H168 configuration, yielding a typical synthesised beam of 
7.8" X 6.3" at half-power beam-width. 

Bandpass calibration was performed on PKS 0537-441 
at the start of every track, and nearby phase calibra- 
tor PKS 0953-1-254 was observed at 10-15 minute intervals 
throughout the observations. The same source was used to 
check the array pointing solution once an hour. PKS 1934- 
638 (the standard flux calibrator at the ATCA) was observed 
at the end of every track to provide primary flux calibration. 

Observations were scheduled at night in order to ob- 
tain the most stable possible conditions for 7mm band ob- 
servations. Conditions were generally good throughout the 
observations, but data taken at elevations < 30deg were 
discarded due to a combination of increased atmospheric 
opacity and less stable conditions at the start and end of 
each night. Images were generated using natural weighting. 
The measured RMS noise at the phase centre in the flnal 
images was typically 74/iJy beam^^ at 37 GHz and 7S fiJy 
beam at 39 GHz in 2MHz (16kms"^) channels. 

A multi-frequency (4GHz bandwidth) synthesis image 
centred at 37.887 GHz was also generated from the data, 
reaching a sensitivity limit (RMS) of 3.1 nJy beam~^. 



3 CONSTRAINTS ON THE Z ~ 8 GRB HOST 

The position of the 2 ~ 8 GRB 090423 is known to sub- 
arcsecond accuracy (from near-infrared imaging of the burst 
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afterglow). The stellar light from galaxies at these dis- 
tances is believed to have a typical half-light radius of 
■ ^O-Skpc (0.06"), b ased on photometrically selected samples 
l|Oesch et al.ll2010l) . The beam size of our ATCA observa- 
tions is two orders of magnitude greater than this, and any 
emission from the the GRB host galaxy itself is likely to be 
contained within a synthesised beam centred at the burst 
location (and probing a region 40 x 30kpc in extent). Our 
initial analysis focused on the properties of our data at this 
location. 

No continuum source was seen at the GRB location in 
the multi-frequency synthesis image, allowing us to place an 
upper limit on the 37.9 GHz continuum flux of this source of 
•S'8mm=9.3 /iJy (3 a, assuming the source is unresolved). This 
observed frequency corresponds to the rest frame 850 pm 
band and so prob es the warm dust in the system. We use 
the prescription of lAravena et al.1 l|2008l ) to calculate an up- 
per limit on the dust mass at z = 8.2, assuming the dust 
emission can be modeled by a modified black body and ac- 
counting for the heating effect of the cosmic microwave back- 
ground at this redshift. A number of significant assumptions 
are required to make this calculation. Most crucially, the 
dust coefficient and emissivity index are fixed at local val- 
ues, which may not be appropriate for high redshift sources, 
particularly young and unevolved sources which are likely 
to be dominated by supernova-processed dust grains (see 
iMaiohno et"al] 12004 IStratta et aIll2007^ . There is also a 
weak dependence on the assumed size of the source. In the 
absence of further data, we assume a source extent of 1 kpc, 
while noting that increasing this by an order of magnitude 
has negligible effect on the derived mass. 

Our derived upper limit on the dust mass in this source 
is 8.0 X 10* Mq at 40K. Given the typical stellar masses of 
GRB host galaxies at low redshift (see section [5]) our non- 
detection of the source in dust continuum is unsurprising. 
The m ean 850 fj,m (observed) flux measured bv lTanvir et al.l 
l|2004h for a sample of 21 GRB hosts at 0.45 < z < 3.42 im- 
plies an average dust mass of only ~ 5 x lO'^ Mq for the lower 
redshift GRB host population using the same prescription 
applied here. For thermal dust emission across a reasonable 
range of temperatures, and for a typical low redshift GRB 
host galaxy, our constraint implies only that the dust mass 
does not exceed the stellar mass. 

In addition to probing the dust continuum of the source, 
our observations are sensitive to CO (3-2) line emission span- 
ning the redshift range 7.7<2<8.6 - easily encompass- 
ing the Az = 0.13 uncertainty in the optical redshift esti- 
mate for the GRB host. We have examined our datacube 
carefully for such emission. Low luminosity galaxies at high 
redshift have been observed to have lower line widths than 
most submillimeter galaxies. Indeed, the only two known 
2 ~ 5, non-AGN, CO-emitting galaxie s at have line widths 
of ~ 160kms~^ (ICoppin et al.1 [2OI0I) and ~ 110 km s"^ 
l|Stanwav et al.ll2008l ) which mav suggest that the molecular 
gas in these sources does not form a rotationally-supported 
disk. With these narrow lines in mind as a model, we select 
150 km s~^ as a plausible velocity width for any emission line 
in our distant source. 

In order to optimise our sensitivity to emission with this 
and a variety of both broader and narrower velocity widths, 
the spectral datacube was rebinned at resolutions of 2, 4, 8 
and 16 MHz and each rebinned cube was analysed indepen- 
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Figure 1. The spectrum of our emission line source (above, see 
section r4.2l l. compared with the non-detection spectrum extracted 
at the location of the GRB host galaxy (below, section [Sjl and 
assuming the latter is a point source. Data are binned to 4 MHz 
(32kms~^) resolution. The shaded region indicates the expected 
frequency range for CO(3-2) emission from the GRB 090423 host 
galaxy, assuming a la error range on the optical/near-infrared 
redshift. The spectra shown span 10,000 km s~^ in velocity space 
and represent only 35% of our total observed frequency range. 
The effect of primary beam attenuation on channel noise at the 
line emitter location is clearly visible. 



dently. Despite this, a spectrum taken through the data cube 
at the location of the GRB host galaxy shows no evidence for 
line emission at the sensitivity limit of the data at any spec- 
tral resolution (see figure [!}. This places an upper limit on 
the velocity integrated line width of the CO (3-2) transition 
(expected at 37.45t^;24 GHz) of ScoAv <0.026 Jy kms^^ 
(3(t) for a 150kms~^ wide emission line. 

The estimation of molecular gas mass from CO line 
intensity requires the application of an empirical calibra- 
tion for Mgas/Lco- For ease of comparison with other au- 
thors, we assume that the molecular gas in this source is 
optically thick and, apply the standard conversion factor 
Mgas/-L'co=0.8MQ (K kms-^)-\ derived for ULIRGS and 
SMGs at intermediate redshifts and believed to be con- 
stant for such sources betwee n 2 = and 2 = 4 (see 
ISolomon fc Vanden Boutll2005l ). Using this conversion we 
determine an upper limit on the molecular gas at z — 8.23 
of M(H2)< 4.3 x 10^ Mq (3 a). 

We note that the conversion has not yet been demon- 
strated for low mass galaxies forming stars at a relatively 
low rate, at significantly sub-solar metallicities metallicities, 
or at early times, for which a higher conversion factor may 
be more appropriate. Applying the Galactic conversion fac- 
tor, Mgaa/ico=4.6 Mq (K kms~^)~^, would yield limits a 
factor of 5.75 times weaker than those presented here. The 
true value of the conversion factor is unknown but likely 
to lie some where between these extremes (see discussion in 
iTacconi et al.„2008i) . 
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4 RADIO SOURCES IN THE GRB FIELD 
4.1 Continuum Sources 

We quantify tlie number density of continuum sources in our 
field, considering a circular region extending from the point- 
ing centre to the radius at which primary beam attenuation 
limits sensitivity to 40% of its maximum value (47"), yield- 
ing a survey area of 5.4 x 10~*deg~'^. We detect a single 
unresolved 5 a source in our 4 GHz-bandwidth continuum 
image, with a peak flux Sgmm = 15.9 ± 3.2 fijy, determined 
by fitting an Gaussian point source to the image data. The 
source lies 33 arcseconds to the south-east of the GRB host 
galaxy. Two other features in the continuum map occur at 
3 a relative to the local noise level. 

In each case, several faint galaxies are visible at the 
source location in available deep J and H band imaging of 
the field (see section 14. 2|) but the large beam size of the 
ATCA at 39 GHz prevents secure identification of the cor- 
rect near-infrared counterpart (if any) or determination of 
an accurate redshift. None of the three continuum sources 
show evidence for strong line emission in the frequency range 
observed. 



4.2 A Molecular Line Emitter 

In principle we are sensitive to any emission line source ap- 
pearing in the telescope primary beam and in our 36-40 GHz 
frequency range. In practice, this region of the spectrum is 
largely devoid of significant line s from Galactic or low red- 
shift sources (|Lovas et al.1119791 ). The only strong features 
seen from Galactic sources in this frequency range (i.e. emis- 
sion lines without comparable-strength neighbouring lines) 
are high order transitions in methanol masers, which are 
usually highly polarised and associated with high mass star 
formation regions. 

There are only three credible extragalactic interpreta- 
tions of an emission line occurring at these high frequen- 
cies. We are sensitive to CO (1-0) line emission occurring at 
z = 2.05 ± 0.15, C0(2-l) at 2 = 5.10 ± 0.30 and CO(3-2) 
line emission at z = 8.15 ± 0.45. 

The data cube obtained contains spectral information 
for any source lying within the primary beam of the ATCA. 
We rigorously inspected our data cube for emission line can- 
didates lying within 47" of the pointing (and hence phase) 
centre. As in the case of the GRB host galaxy, we consider 
the data rebinned in frequency on scales ranging from 16 to 
126kms~^ per spectral resolution element, searching both 
for particularly narrow line emission and for weaker but 
broader emission lines. Clear detections at the same spatial 
location in three or more adjacent channels were required to 
identify an emission line candidate. While pattern noise in 
the reduced images often creates strong peaks in the noise 
distribution, these are limited to one spectral channel, and 
noise peaks in adjacent channels appear spatially offset from 
one another. 

We identify a single strong emission line source in 
our spectral data cube, centered at 37.164 GHz as shown 
in figure [1] The best fitting Gaussian line profile has a 
peak fiux of 0.58 ± 0.12 mjy and a velocity full width at 
half-maximum of 96kms~^, yielding an integrated line fiux 
Sco^v = 0.059 ± 0.012 Jykms"^, accounting for primary 



beam attenuatiorlf|. The line flux is not signiflcantly po- 
larised (although we note that at this low signal to noise 
there is large fiux uncertainty in each polarisation). The 
source, located at 09''55"'35=96 +18°08'24"7, is not detected 
in 8mm continuum emission. 

As discussed above, three plausible redshifts exist for 
this source: z — 2.10, 5.20 and 8.31, corresponding to the 
first three rotational transitions of the carbon monoxide 
molecule. Due to the offset between redshift, luminosity dis- 
tance and intrinsic line flux for these three transitions, our 
observations reach comparable mass sensitivities in all three 
lines: 1.0, 1.1 and 0.98 xlO^'^M© in molecular gas respec- 
tively (assuming Mga^/ L'qq—0.8 as discussed in section|31). 

Since the galaxy almost certainly does not flU our 7.8" x 
6.3" beam, a detailed calculation of the dynamical mass is 
not possible with the current data, and, in fact, it is not 
clear whether the line width is actually diagnostic of ordered 
rotation in sources of this type. However, while narrow, the 
emission line velocity width is consistent with the calculated 
gas mass if the molecular gas is distributed in a virialised 
disk one arcsecond in diameter at 2; = 8 (0.6 arcseconds at 
= 2), seen close to face-on. 

Given the hierarchical model for galaxy formation, and 
in the absence of additional information on large scale struc- 
ture in this field, our a priori most likely interpretation has 
to be that this is an intermediate-mass star forming galaxy 
at z = 2.10. Such a galaxy would typically be easily de- 
tec ted in deep optica l (7 ~ 24 in similar mass examples, 
see ICasev et al.ll2009l ) and near-infrared data, although the 
presence of significant amounts of dust could lead to the 
absence of a counterpart (as is the case for submillimetre 
luminous galaxies). Higher redshift sources would typically 
be both fainter and bluer in the near-infrared. 

In figure [2] we show existing deep near-infrared imag- 
ing of the region from HST/WFC3 (P.I. Tanvir, analysis in 
prep). This data reaches a limit of J ~ 28.5 (AB, 3 cr) in a 1" 
aperture and was obtained in order to study the host galaxy 
of GRB 090423. Near-infrared data of this depth and qual- 
ity is rare. Figure [5] demonstrates the high-surface density of 
faint sources in such data, and hence the difficulty of assign- 
ing an unambiguous counterpart to the line emitter. The am- 
biguity is further heightened by uncertainties in the pointing 
accuracy of the ATCA. At millimetre wavelengths, the tele- 
scope is used in reference pointing mode, with a nearby cali- 
brator (PKS 0953-1-254 in this case) used to check and refine 
the pointing solution once an hour throughout observations. 
Typical corrections at the end of each hour's observations 
were of order a few arcseconds, and rarely exceeded half the 
size of the synthesized beam. 

Given this uncertainty, there are five near-infrared 
sources, well-detected in both bands, which may represent 
a counterpart for the line emission source, as well as a 
number of low-surface brightness features at the noise level 
of the data. The J — H colours of these sources are in- 
sufficient to provide a redshift estimate, although we note 
with interest that the faint unresolved point source clos- 
est to the centre of the line emission (marked '1' in fig- 
ure [2]) shows very blue near-infrared colours (J — 27.4, 



^ which reduces the measured flux at this position to 42.5% of 
its true value 
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Figure 2. Deep near-infrared imaging in the F125W (J) band and F160W (H) band, overlayed with contours showing the spatial position 
of CO line emission, integrated across the 22 MHz width of the line emission. Contours mark 5, 7.5 and 10 cr levels relative to the noise 
in the line-integrated image. Images are 12 arcseconds on a side and oriented with north uppermost. As discussed in section 14.21 the 
synthesised beam of our observations at this frequency is 8''00 X 6'.'45 and the pointing is accurate to approximately half a beam (roughly 
corresponding to the second contour). It is apparent that a number of plausible counterparts exist for the line emitter with the five 
best-detected labelled in the J band image 



J — H = —0.83 ± .15), a trait typical of hiRh redshift 
(z > 5) galaxies (e.g . IStanwav. McMahon. fc Bunked l2005l : 
iBouwens et al.|[2010l i and difRcult to obtain for normal stel- 
lar populations at lower redshift. 

Further data will be needed to unambiguously confirm 
the redshift of this source, and in the absence of such data, a 
tentative identification as CO (1-0) emission from a faint z = 
2.10 starburst galaxy remains the most likely interpretation. 



5 DISCUSSION 

5.1 Constraints on tlie GRB Host Galaxy 

A typical L* Lyman break galaxy at g ~ 7 — 8 is believed to 
have a stellar mass of a few x 10^ Mq (|Gonzalez et al.ll2010l : 
iLabbe et al.ll2O10l ). although considerable uncertainties (pri- 
marily stellar population synthesis model used, metallicity, 
photometric uncertainty and the deblending of faint, con- 
fused sources) remain in this calculation. This is very com- 
parable to the median stellar mass de termined for GRB hos t 
galaxies at z < 1.2 (1.3 x 10^ M©, ISvensson et al.ll2010l) . 
and little evolu tion in this is seen to much higher redsh ifts 
either directly jSavaglio. Glazebrook. &: LeBorgnell2009l ) or 
in semi-analytic simulations tuned to match the observed 
distribution at z = 1 — 6 (Chisari, Tisscra, & PcUizza 20ic|). 

Our constraint on the mass of molecular gas in the host 
galaxy of GRB 090423 requires Mh2 to be within a factor of 
a few of this expected stellar mass for GRB hosts, suggesting 
either that the GRB host is less massive than those typical 
at lower redshift or that its star formation is already well 
advanced, i.e. that it must be at least a third of the way 
through its starburst process, assuming 100% efficiency in 
the conversion of gas to stars. 

This tight limit is sensitive to gas not only in the host 
galaxy itself, but in the 40 x 30kpc region encompassed by 



the synthesised beam of the observations. As a result the up- 
per limit includes gas infalling on or recently ejected from 
a star forming galaxy, and in any close neighbours. As such 
it is somewhat surprising that a galaxy that is likely the 
progenitor of a much more massive low redshift system has 
such a limited fuel supply available for later and ongoing 
star formation. While no data is currently available on the 
current star formation rate of this source, it is instructive to 
consider a crude model in which the host galaxy is typical 
of lower redshift GRB hosts (i.e. M* Ri 10^ M©) and star 
formation began at 2; ~ 12 (during the epoch of reionisa- 
tion), yielding a a stellar population age of 240 Myr and a 
continuous star formation rate of 4M0 yr~^ (comparable to 
th e median of 3.8M(;)y r~^ found in the GRB host sample 
of lSvensson et al]l2010^ . In this model, the host galaxy will 
exhaust its fuel supply hy z ~ 5.7 at the current rate of star 
formation. If the host is less massive or star formation began 
earlier (and hence the average star formation rate is lower) 
the star formation epoch may extend to lower redshifts. 

We note that this conclusion is drawn from the as- 
sumption of an optically-thick thermalised, metal-enriched 
gas, and from a CO to H2-|-He gas mass conversion factor 
that is calibrated on the most massive starbursts at inter- 
mediate redshifts. These assumptions may well be flawed 
on a number of levels. As mentioned in section |3l the ap- 
plicability of a simple conversion (widely known as the X- 
factor) from CO (1-0) luminosity to a gas mass has not been 
demonstrated at higher redshifts, or in galaxies with sub- 
ULIRG star formation rates. The low Mgas/Lco ratio in 
ULIRGs is believed to arise from the presence of a galaxy- 
wide starburst in which the bulk of the interst ellar medium is 
both dense and irradiated by star formation l|Solomon et al.l 
Il997h . If the GRB host galaxy at z ~ 8 is, like those at 
intermediate times, similar in character to 2 5 Lyman 
break galaxies in terms of mass, metallicity and star forma- 
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tion rate, then it is interesting to note that such systems do 
indeed seem to be gal axy-wide s tarbursts with high specific 
star formation rates (jVerma e t al- ,2007. ) and no evidence 
that the UV-luminous knots observed are merely super- 
starclusters embedded in a more mass ive, cooler medium 
ijPavies et al.ll20ld : [Stanwav et al1l2010l l. While this may ar- 
gue in favour of the ULIRG-like conversion ratio used in this 
paper, the low metallicity and comparitively low star forma- 
tion rates of such sources may have the effect of increasing 
the ratio. This significant source of uncertainty must neces- 
sarily be addressed by later studies. 

Although GRB 090423 is believed to be a population II 
event and thus indicative of a metal enriched star forming 
population, the fraction of relatively pristine II2 gas in the 
outer regions of a galaxy at early times is likely to be higher, 
particularly well outside the star forming disk. Our observa- 
tions cannot rule out the presence of a large reservoir of such 
gas, either in the host galaxy or in filaments surrounding it 
since this would likely not be reflected in the measured CO 
luminosity. 

Similarly, the conversion assumes a known ratio of emis- 
sion in the higher order rotational lines to that in the CO(l- 
0) fundamental transition. At high redshifts, the cosmic mi- 
crowave background provides a minimum temperature for 
a typical galaxy, pushing the rotational excitation ladder 
of carbon monoxide to higher levels and hence reducing 
the measured flux in the lower excitati on states. The ef- 
fect of this is seen in the recent models of lObreschkow et al.l 
l|2009bl ). which indicate that the highest surface density of 
sources to a given flux limit is seen in CO (3-2) emission 
at low redshift and in CO (5-4) emission by 2 = 5. Future 
follow-up of this and similar sources with ALMA (which, 
while having a more limited fleld of view, will access higher 
order CO transitions) may yield an indication of this effect 
and prove more successful in determining the precise source 
redshift and gas mass. 



5.2 Large Scale Structure at High Redshifts 

As mentioned in section our deep datacube probes cos- 
mological volume in each of three redshift ranges. Consider- 
ing the region in which our flux sensitivity is >40% of that at 
the pointing centre yields a field of view of 5.4 x 10"* deg"'^, 
while the 4 GHz bandwidth yields a line of sight depth of 
Az = 0.29, 0.59 and 0.88 at z = 2, 5 and 8 respectively. 

Given a single emission line source (with redshift yet 
to be determined) we are able to place an upper limit on 
the surface density of CO-luminous galaxies in each of these 
redshift ranges. Our constraints are shown in figure |31 to- 
gether with the predicted surface densities of sources with 
L'qq > 0.05 Jykms"^ (i.e. to a fiux limit matching our de- 
tection of a sin gle possible line) de t ermine d from the S-cubed 
simulations of lObreschkow et al.l l|2009al ). These analyse a 
conical section through the Millennium Simulation, applying 
a semi-analytic prescription to calculate the predicted CO 
rotation emission line excitation ladder and H2 gas proper- 
ties associated with about 2.8 x 10* dark matter haloeqj. As 
figure 13] demonstrates, very few z = S galaxies in the simula- 
tion would satisfy our fiux selection limit in the CO (3-2) line 
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Figure 3. Our upper limit on the number density of line emis- 
sion sources satisfying L'^q > 0.05 Jykms"^, given the velocity 
integrated flux of our single detected line emitter (i.e. assuming 
that we detect 1 or fewer sources in each redshift range). Also 
shown are predicted number densities of sources satisfying the 
same constraint in each emission line from the semi-analytic mod- 
els of lObreschkow et al. I l l2009al') as discussed in section [5]2l Since 
these models are based on a section through the millennium simu- 
lation, they are subject to small number statistics in the predicted 
number densities, particularly at the highest redshifts. 



(the simulation field of view at z = 8 is about 4x4 degrees^), 
and hence the uncertainty on the predicted number densities 
is dominated by small number statistics. 

Interestingly, regardless of which line identification is 
considered most likely, our detection of a CO-emitting source 
is something of a surprise. Given these simulations, and in 
the absence of lensing, we might expect to have detected 
either a, z — 2.1 CO(l-O) emission line or a 2; = 5.2 CO(2- 
1) emission line of this strength in one of 35 blind ATCA 
pointings. It is possible that one of the intermediate redshift 
galaxies visible in figure [2] is a lensing source, boosting the 
flux ofa.z = 2oTZ — 5 background galaxy, or that we have 
chanced across rare large scale structure at these redshifts. 

The probability of a randomly selected pointing detect- 
ing a, z — 8 galaxy in CO emission is considerably lower 
(approximately one in two thousand). However this fleld is 
not randomly selected. In the unlikely event that follow-up 
investigations confirm the very blue {J — H — —0.8) faint 
galaxy as a 2 = 8 near-infrared counterpart to the emission 
line source, this field would host both a Gamma Ray Burst 
and a more massive gas-rich source (separated by a pro- 
jected distance of 200 kpc), representing a highly significant 
overdensity above the typical galaxy distribution at this red- 
shift and suggesting that GRBs are indeed effective tracers 
of large scale structure at high redshift. While constraints 
from a single GRB and associated field remain weak, such a 
detection might argue against the suppression of star forma- 
tion in low m ass haloes seen in warm dark matter models, as 
discussed by iMesinger. Perna. fc HaimanI 1.2005') which are 
expected to yield both fewer high redshift GRBs and fewer 
star forming galaxies associated with them. 



Even by 



8, the effects of cosmic variance can be 



http: / /s-cubed. physics. ox. ac.uk/s3_sax/sky 



significant, (ihev et all l|2010l ) have suggested a two order 
of magnitude variation in the density of 0.5 Mpc^ regions 
at this redshift, based on large N-body simulations probing 
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structure formation during reionisation. Given the signifi- 
cant expenditure of telescope time necessary to obtain deep 
imaging suitable for exploring the z ^ 8 Universe, exploring 
underdense regions is likely to be both disappointing and 
uninformative. Following up GRB fields (should additional 
z > 6 bursts be identified) may well be the most efficient 
method for identifying the high density regions that are the 
sites for later massive galaxy and cluster formation. 



6 CONCLUSIONS 

We have observed the host galaxy of GRB 090423 and its 
immediate environs at radio frequencies (36.1-39.7 GHz), in 
order to search for molecular gas emission in the CO (3-2) 
line. Our main conclusions can be summarised as follows: 

(i) No continuum source was seen at the GRB loca- 
tion in the multi-frequency synthesis image, allowing us to 
place an upper limit on the 37.9 GHz flux of this source of 
S'8mm=9.3 ^Jy {3a, assuming the source is unresolved). 

(ii) No line emission is detected from the GRB 090423 
host galaxy. Applying standard conversion factors, we de- 
termine an upper limit on the molecular gas in this z — 8.23 
source of M(H2)< 4.3 x 10^ M© {3 a). 

(iii) This tight limit on the fuel for star formation sug- 
gests either that the GRB host is less massive than those 
seen at lower redshift, or that its starburst is well advanced. 
We caution that this conclusion depends critically on the 
assumptions required to calculate gas masses, which are not 
yet well calibrated at high redshift. 

(iv) We identify a single strong emission line source in 
our data, centered at 37.164 GHz and 44 arcseconds to the 
south-west of the GRB. The emission has a peak flux of 
0.58 ± 0.12 mjy and a velocity full width at half- maximum 
of 96kms~^. This corresponds to a molecular gas mass of 
1.0, 1.1 or 0.98 xlO^° Mq if interpreting the line as CO(l-O) 
at z = 2.1, CO(2-l) at z = 5.2 or CO(3-2) a,t z = 8.3 respec- 
tively. Efforts to identify and characterise a counterpart to 
this emission line source are continuing. 



NOTE ADDED IN PRESS 

Flux limits given in sections 2 and 3 are obtained from data 
that has been 'cleaned' using MIRIAD. Further inspection 
and consultation with ATCA staff suggests a less vigorous 
cleaning strategy may be appropriate. This would result in 
an effective increase in the image noise, and hence quoted 
flux and mass constraints, by a factor of two (~3 if clean- 
ing is omitted completely from the data reduction). Quoted 
fluxes for the line emission candidate are unaffected, and 
there is no substantive impact on our discussion or science 
conclusions. 
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